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SECTION  I 


INTRODUCTION 

This  report  presents  the  results  of  an  investigation  to  develop  an 
analytical  method  which  predicts  the  maximum  pressure  distortion  level 
and  provides  a  synthesized  pressure  distortion  map  at  the  entrance  to  the 
turbine  engine  compression  system.  The  method  uses  the  steady-state 
total  pressures  and  the  statistical  properties  of  the  time-variant  total 
pressures  measured  at  the  compressor  face  plane. 

A  high  degree  of  compatibility  between  the  inlet  and  engine  is 
necessary  for  an  aerodynamically  stable  propulsion  system.  An  essential 
element  in  assessing  that  compatibility  is  the  determination  of  the 
tinic-variant  total  pressure  distortion  generated  by  the  inlet.  A 
deterministic  analysis  procedure,  shown  in  Figure  1,  has  evolved  for 
processing  time-variant  pressure  distortion  data.  Filtered  high-response 
compressor  face  pressure  data  are  combined  with  their  respective  steady- 
state  total  pressure  component  and  input  to  a  set  of  engine  distortion 
parameter  equations.  Peak  distortion  levels  and  pressure  distortion  maps 
are  determined.  If  the  peak  distortion  level  is  less  than  some  limit 
value  related  to  an  engine  surge  margin  allocation,  a  compatible 
inlet-engitie  combination  exists.  The  procedure  has  been  a  successful  but 
conservative  approach  to  the  problem. 

Analog  editing  systems  have  been  used  extensively  to  provide  that 
assessment  of  inlet-engine  compatibility.  One  such  system  called  DYNADFC 
(Dynamic  ^ata  Siting  and  Computing),  utilizes  a  modern  state-of-the-art 
hybrid  computer  (Reference  1).  DYNADEC  provides  for  the  continuous 
solution  of  the  distortion  parameter  equations  and  the  handling  of  the 
massive  quantities  of  data  in  a  timely  and  cost-effective  manner.  Entire 
data  records  can  be  screened  in  a  real-time  environment  with  final 
results  available  immediately.  The  high  degree  of  system  flexibility 
allows  program  parameters  to  be  easily  changed  and  thus  provides  for  an 
efficient  analysis  process. 

While  the  analysis  process  appears  to  be  well  in-hand,  the  acquisi¬ 
tion  of  dyriamir  pres'^ure  data  has  become  an  increasingly  expensive 
proposition  ciue  to  higher  wind  tunnel  operational  costs,  the  use  of  large 


To  illustrate  the  limitations  ct  the  basic  model,  the  cifference 
between  the  predicted  and  measured  peak  distortion  level  as  a  function  of 
measured  distortion  level  is  shown.  Next,  the  average  starioard  deviation 
of  the  difference  between  predicted  and  measured  probe  pressure  recoveries 
as  a  function  of  turbulence  is  presented.  The  °f'PKSp,pp  for  a  single  case 
is  defined  as  o  i 


n  -  1 


E—  -  — 


where  aPj  is  the  difference  between  a  predicted  and  measured  pressure  and 
aPj  is  the  mean  difference  between  the  predicted  and  measured  pressures. 
Finally,  the  predicted  and  measured  pressure  contour  maps  are  shown, 
including  the  distortion  levels  and  the  histogram  of  the  difference 
between  the  predicted  and  measured  probe  pressure  recoveries. 

Figures  4  through  8  present  predicted  and  measured  results  for  the 
GE^  engine  distortion  parameters.  The  measured  peak  distortion  levels 
and  pressure  distortion  maps  were  determined  from  screened  inlet  pressure 
data  filtered  at  a  cut-off  frequency  of  500  HZ.  The  differences  between 
predicted  and  pleasured  peak  distortion  levels  are  shown  in  Figure  4. 
Peak  distortion  levels  are  expressed  in  terms  of  APRSp  with  the  circum¬ 
ferential  and  radial  components  being  and  respectively. 

The  basic  model  significantly  overpredicts  APRSp  and  The  radial 

component  is  more  randomly  distributed,  but  with  some  large  differences 
as  well. 

In  Figure  5,  the  average  standard  deviation  of  the  difference  ir 
probe  pressure  recoveries  is  presented.  The  difference  between,  predicted 
and  measured  pressures  increases  with  increasing  average  turbulence. 
Since  '.iRMSpjpp  will  vary  with  different  sets  ot  random  numbers,  several 
cases  were  repeated  several  times  using  different  starting  seed  values. 
A  litiear  regression  analysis  (Reference  15)  was  then  accomplished  on  the 
total  set  of  account  tor  differences  due  the  raridom  number  process.  'he 
resulting  regression  line  and  others  that  foi'iow  serve  as  baselines  Ic 
show  improvements  to  the  model. 

Figures  6,  7,  arc  &  compare  the  predicted  and  measured  peak 

distortion  maps  lor  three  levels  of  average  turbulence.  For  all  three 


1  fi 


defined  as  they  are  not  a  part  of  the  input  to  the  statistical  prediction 
method. 

Compressor  face  instrumentation  used  to  measure  the  time-variant 
pressures  consisted  of  forty  steady-state  and  high  response  probes  in  an 
eight-rake  by  five-ring  array.  The  fluctuating  pressure  data  was  filtered 
at  500  HZ  and  1000  HZ  (-3db)  consistent  with  engine  sensitivity,  inlet 
model  scale,  and  available  cut-off  filter  frequencies. 

The  Pratt  &  Whitney  Aircraft  (PWA)  K,  and  the  General  Flectric 
APRSp  (GEj)  and  IDL  (GE2)  were  the  three  sets  of  engine  distortion 
parameters  used  in  llie  analysis.  For  the  PWA  distortion  methodology,  the 

K 

total  distortior.  is  the  sum  of  the  ci rcumferential  distortiori,  K., 

^  K  ^ 

iHK*  the  weighted  radial,  distortion,  ^RA... 

C 

\  ^  ,2, 

The  GE  parameters  describe  either  the  stall  pressure  ratio  loss,  APRSp, 
the  sum  of  the  loss  in  fan  stall  pressure  ratio  due  to  circumferential 
and  redial  distortion,  or  the  ratio  of  surge  margin  required  to  that 
available  for  distortion  defined  as  IDL.  Functionally,  APRSp,  and  IDL 
can  be  expressed  in  terms  of  the  maximum  circumferential  and  radial 
distortion, 

APRSp,  IDL  -  f  (IDC^^^,  IDP^^^)  (3) 

Basic  Nodel  Results 

A  common  format  of  comparison  between  preoicted  and  measured  quanti¬ 
ties  is  presented  throughout  this  section  and  in  Section  IV.  Several 
measures  of  goodness  are  used  to  assess  the  capability  of  the  basic  end 
improved  models.  These  measures  of  goodness  were:  the  difference  between 
the  predicted  and  measured  maximum  distortion  level,  the  average  standard 
deviation,  1  he  distribution,  and  the  range  of  the  difference  between 
predicted  and  measured  pressures,  the  similarity  between  the  predicted 
oressure  and  measured  contour  maps,  and  tlie  number  of  predicted  pressures 
within  .1  percent  of  their  measured  values.  Peak  distortion  pressure 
contour  maps  determined  with  DYNADEC  serve  as  the  basis  for  the  measured 
va 1 ues . 
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Figure  3  Basic  Statistical  Prediction  rV)del 


SECTION  III 


BASIC  MODEL  ANALYSIS 

Statistical  Predict ion  Model  Descripti or 

The  basis  for  the  method  is  that  a  synthesized  fluctuating  presbu^e 
component  can  be  constructed  and  added  to  the  steady-state  pressure  to 
form  the  dynamic  total  pressure.  The  fluctuating  pressure  is  assumed  to 
be  stationary,  uncorrolated,  and  random  with  a  ttonnal  distribution. 

The  synthesized  fluctuating  pressure  can  be  determined  for  each 
probe  using  a  random  number  generator  with  a  zero  mean,  and  the  standard 
deviation  derived  from  the  measured  turbulence  levels.  If  the  dynamic 
total  pressure  is  defined  in  terms  of  pressure  recovery,  then  for  each 
probe: 


The  basic  model,  Figure  3,  consist?  of  two  fundamental  elements, 
the  generation  the  compressor  face  dynamic  total  pressures  and  the 
determination  of  the  maximum  distortion  level  and  pressure  contour  map. 
PMS  turbulence  anu  random  numbers  are  combined  to  form  the  fluctuating 
pressure  components  which  are  added  to  the  steady-state  pressures  to  fomi 
the  dynamic  total  pressures.  The  dynamic  total  pressures  are  input  to 
tile  distortion  parameter  equations  and  the  level  is  determined.  The 
pressure  contour  rrap  is  also  generated.  Forty  new  random  numbers  are 
geru-raied  providing  a  new  set  of  dynamic  total  pressures  which  represent? 
delta  from  another  equivalent  time  slice.  The  distortion  level  is  computed 
and  compared  to  the  current  maximum  value.  The  sequence  is  repeated 
crtil  a  desired  sample  size  is  reached. 

JjT_let  Data  Base  and  Engine  Distortion  Parameters 

Forty-nine  cases  from  two  sets  of  inlet  data  with  three  sets  of 
engine  distortion  parameters  were  used  in  the  investigation.  Average 
compressor  face  turbulence  ranged  from  approximately  .01  to  .08.  Inlet 
opera' ii,g  conditions,  test  conditions,  and  inlet  gec'tretry  have  not  been 


A  third  investigator  to  report  on  the  development  of  a  statistical 
synthesis  method  using  a  random  number  generator  and  local  RMS  pressure 
data  was  Borg  (Reference  2).  In  his  method,  the  normally  distributed 
random  numbers  were  generated  by  adding  twelve  independent  random  numbers 
from  a  uniform  distribution  using  the  central  limit  theorem.  Briefly, 
the  central  limit  theorem  states  that  the  sums  of  independent  random 
variables  under  fairly  general  conditions  will  be  approximately  normally 
distributed,  regardless  of  the  underlying  conditions  (Reference  14).  The 
RMS  pressures  were  determined  using  a  suitable  cut-off  filter  frequency. 
Borg  used  the  Pratt  &  Whitney  and  Rolls  Royce  distortion  parameters  in 
evaluating  his  method.  Results  with  the  PWA  parameters  were  consistent 
with  that  of  Stevens  and  Sanders.  One  interesting  aspect  of  the  study 
was  the  use  of  average  RMS  pressure  values  at  all  probe  locations.  Since 
there  was  no  drastic  change  in  the  distortion  level  correlations,  Borg 
concluded  that  a  reduced  number  of  turbulence  values  could  be  used  if 
they  reflected  the  average  turbulence  intensity.  No  comparisons  were 
made  between  predicted  and  measured  peak  distortion  maps. 

In  summary,  the  statistical  prediction  methods  were  found  to  provide 
0  reasonable  prediction  of  the  maximum  distortion  level.  The  methods 
developed  by  Motycka,  Stevens,  and  Melick  also  provided  a  synthesized 
pressure  distortion  map.  Each  of  these  methods  were  found  to  predict  a 
reasonable  distortion  map  at  low-turbulence  levels,  however,  as  the 
turbulence  level  increased  and  as  the  number  of  measured  values  of 
turbulence  used  in  the  analysis  decreased,  the  quality  of  the  pressure 
distortion  map  deteriorated. 

These  investigators  proposed  the  use  of  statistical  methods  to 
achieve  cost  reductions  for  the  analysis  of  time-variant  distortion  data 
from  small  models  with  limited  instrumentation.  Such  methods  were  not 
intended  to  supplant  the  editing  systems  used  for  screening  inlet 
distortion  data  during  validation  tests  of  inlet-engine  compatibility  of 
aircraft  systems. 
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expected  engine  sensitivity  in  determining  the  probe  RMS  turbulence 
levels.  Second,  no  power  spectrum  shaping  is  used  in  the  analysis. 
Large  sample  sizes  are  considered  necessary  to  accurately  represent 
dynamic  pressure  data.  For  example,  if  the  highest  frequency  of  interest 
is  approximately  1000  HZ,  and  five  or  more  samples  are  required  per 
cycle,  then  some  150,000  samples  or  time  slices  are  required  to  analyze 
an  equivalent  30  seconds  of  data. 

Predicted  peak  distortion  levels  and  pressure  contour  maps  were 
compared  to  an  F-15  and  F-18  inlet  data  base.  Peak  distortion  levels, 
based  on  Pratt  &  Whitney  and  General  Electric  distortion  parameters,  were 
approximately  10  percent  under  and  over  the  measured  values,  respectively. 
This  result  is  consistent  with  Sander's  analysis  of  Motycka's  method. 
Stevens  also  investigated  the  use  of  fewer  measured  turbulence  values  to 
predict  the  peak  distortion  levels  and  found  that  the  results  were  about 
the  same  as  when  all  the  measured  RMS  turbulence  values  were  used.  A 
comparison  of  predicted  versus  measured  pressure  contour  maps  using  both 
sixteen  and  forty-eight  RMS  turbulence  values  was  also  accomplished.  It 
was  concluded  that  the  pressure  contours  obtained  using  the  statistical 
model  with  forty-eight  RMS  turbulence  values  generally  agreed  well  with 
the  measured  contours,  but  for  the  reduced  number  of  turbulence  values, 
the  predicted  contours  were  only  representative  of  the  measured  contours. 

Forner  and  Manter  (Reference  13)  used  Stevens'  method  to  predict 
peak  distortion  levels,  based  on  the  Williams  Research  Corporation 
distortion  parameters,  for  cruise  missile  inlet  configurations.  Compar¬ 
isons  between  predicted  versus  measured  peak  distortion  levels  showed 
that  almost  all  of  the  data  fell  within  a  +10  to  -20  percent  band.  For  a 
number  of  cases  there  was  poor  agreement  between  the  predicted  and 
measured  pressure  distortion  maps.  Other  patterns  were  judged  to  be  good 
on  a  qualitative  basis,  but  exhibited  fairly  large  differences  between 
individual  predicted  and  measured  probe  values.  The  duct  flow  was  known 
to  be  highly  dynamic  with  regions  of  separated  flow.  Forner  and  Manter 
olso  investigated  the  use  of  fewer  turbulence  measurements  and  lound  that 
four,  eight,  or  twelve  values  proviaed  as  good  a  prediction  of  distortion 
level  as  with  forty  turbulence  measurements.  As  fewer  measured  turbulence 
values  were  used,  however,  fewer  predicted  patterns  agreed  with  measured 
patterns . 
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Motycka's  method  (Reference  3)  consists  of  determining  the  mean  RMS 
pressure  as  a  function  of  frequency,  and  the  power  spectral  density  for 
each  pressure  probe.  An  amplitude  probability  density  curve  is  generated 
for  each  probe  from  the  RMS  level  assuming  a  normal  distribution.  Random 
numbers  ar(>  converted  to  pressures  from  a  cumulative  amplitude  probabil¬ 
ity  density  function  determined  from  the  integration  of  the  amplitude 
probability  density  curve.  The  synthesized  pressures  are  scaled  to  the 
experimental  PSDs  with  a  digital  filter  and  stored.  Filter  coefficients 
are  determined  from  an  amplitude  gain  curve  formed  by  dividing  the  PSD  of 
the  test  data  by  the  PSD  of  the  random  numbers.  Therefore,  the  power 
spectrum  for  the  synthesized  pressures  are  modified  to  have  the  same 
power  spectrum  of  the  experimental  data.  The  resulting  equivalent 
pressure-time  traces  are  then  reduced  in  the  same  manner  as  digitized 
pressure  data  used  in  the  deterministic  method  for  finding  the  maximum 
distortion  level  and  pressure  contour  map. 

Motycka  examined  a  single  case,  and  while  there  was  good  agreement 
between  measured  and  predicted  results,  he  offered  no  general  conclusions 
regarding  the  accuracy  of  the  method.  Using  the  extreme  value  analysis 
of  Reference  12,  Motycka  concluded  that  the  influence  of  the  filter  was 
to  reshape  the  extreme  value  distortion  versus  tinie  telationship.  As  a 
consequence,  he  suggested  that  the  digital  filter  could  be  eliminated  if 
the  maxiniun  expected  value  is  for  a  relatively  long  inlet  operating  time. 

Sanders  (Reference  10)  evaluated  Motycka 's  method  using  three 
(iistortion  factors  anc  four  sets  of  inlet  data.  There  was  a  tendency  to 
overpredict  the  maximuni  distortion  level  using  General  Electric  (GE) 
distortion  parameters,  and  a  tendency  to  underprediit  the  measured  values 
using  Pratt  8  Whitney  Aircraft  (PWA)  and  Williams  Research  Corporation 
(WRC)  distortion  parameters.  This  tendency  of  underpredicting  was 
believed  to  be  due  +c  an  invalid  assumption  of  normality  for  the  pressure 
data.  The  prediction  of  distortion  patterns  was  considered  good,  but  the 
agreement  between  measured  and  predicted  maps  decreased  signif icantly  for 
turbulence  levels  greater  than  .02. 

Ihe  method  developed  by  Stevens,  Spong,  and  Cliphant  (Reference  4) 
is  very  similar  to  that  of  Motycka,  with  two  differences.  First,  ttie 
pressure  data  is  processed  at  a  cut-off  filter  frequency  consistent  with 


The  results  indicated  that  twenty  probes  yielded  essentially  the  same 
results  as  forty  probes,  with  four  probes  yielding  a  variation  of  less 
than  10  percent  compared  to  a  predicted  distortion  value  based  on  forty 
probes. 

Sanders  (Reference  10)  provided  a  comprehensive  analysis  of  MelitT's 
method  using  three  different  distortion  factor  methodologies  and  four 
sets  of  inlet  data.  One  particular  distortion  parameter  was  modified  by 
using  empirically  derived  statistics  instead  of  theoretically  derived 
statistics.  In  general,  an  excellent  correlation  was  found  for  the  index 
with  the  empirically  derived  statistics,  with  poor  to  good  agreement  for 
the  other  distortion  factors  that  were  based  on  the  theoretical  derived 
statistics.  Comparisons  were  made  of  measured  and  predicted  maximum 
distortion  patterns  for  a  narrow  range  of  low-turbulence  values 
( .01?'T<.0?3) .  Ouo Illative  agreement  appeared  to  range  from  poor  to  good 
for  these  1  ov.'-turbulence  cases.  Schweikhard  (Reference  11)  also  made 
comparisons  n.easured  data  with  Melick's  method  and  found  much  the  same 
results  tor  both  the  distortion  level  and  pressure  distortion  pattern. 

The  next  tiuee  methods  employ  a  random  number  process.  The  funda¬ 
mental  premise  o^  these  methods  is  that  a  random  number  process  can  be 
used  to  synthesi i*'  the  fluctuating  component  of  the  dynamic  total  pressure 
from  the  s tat  i s *■  ica  1  properties  of  the  inlet  pressure  aata.  The  dynamic 
tote'  pressure,  ^T,  consists  v"'  '^wo  components,  the  steady-state 

n  byn ,  p 

iot.ponent,  ,  and  the  fluctuating  component,  T ,  ,  such  that 
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where  T.he  f  luf  ri.atiiu;  ci  opotient  nas  zero  mean  vatue.  It  is  assumed 
that  the  t irne-var icr  f  pressure  data  are  rardom,  stationary,  with  a  normal 
diitnbutioti.  A  randon  tiuiiiber  generator  with  a  normal  distribution  and  a 
zero  inear,  and  the  n.tcsured  standard  deviation  of  pressure  are  used  to 
'’uriii  d  synthesized  fluctuating  pressure  component.  The  dynamic  pressure 
nodel  is  depicted  in  Figure  ?.  The  dynamic  pressures  are  then  input  to  a 
set  of  distortion  parareters  and  a  maximum  value  is  determined. 
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of  a  randoiti  distribution  of  discrete  vortices  being  convected  downstream 
by  the  mean  flow  and  that  the  effects  of  a  series  of  vortices  can  be 
determined  by  summing  the  individual  effects  of  eacn  vortex. 

The  method  assumes  that  most  inlet  total  pressure  fluctuations  are 
normally  distributed  and  that  any  distortion  parameter  computed  from 
these  pressures  will  similarly  be  normal.  The  vortices  in  the  flow  are 
of  arbitrary  size,  strength,  location,  and  orientation.  Vortex  size  and 
strength  are  represented  by  a  Beta  distribution,  while  the  location  and 
orientation  are  assumed  to  be  uniformly  distributed.  Certain  statistical 
functions  are  assumed  to  be  related  to  the  physical  characteristics  of 
the  turbulent  flow  model.  For  example,  the  mean  square  of  the  vortex 
strength  is  proportional  to  the  root-mean-square  of  the  fluctuations,  the 
mean  core  size  is  defined  by  the  frequency  of  the  power  spectral  density 
turiction,  and  the  vortex  flux  is  related  to  the  probability  density 
function. 

Through  a  number  of  assumptions,  the  complex  statistical  and  fluid 
flow  relations  are  simplified  to  a  graphical  solution.  This  graphical 
approach  relates  the  standard  deviation  of  the  distortion  parameters  to 
the  standard  deviation  of  the  pressure  fluctuations,  to  the  mean  vortex 
core  size,  and  to  the  vortex  flux.  The  probability  density  function  of 
each  distortion  factor  must  also  be  defined,  and  a  Beta  density  function 
is  constructed  to  satisfy  constraints  concerning  the  shape  of  the 
distribution. 

To  create  the  instantaneous  distortion  pattern,  a  linear  vortex  is 
located  across  the  duct  between  the  high-  and  low-pressure  regions  of  the 
steady-state  distortion  pattern.  A  total  pressure  increment,  due  to  the 
total  pressure  variation  within  the  vortex,  is  added  to  each  probe 
pressure  value.  The  strength  of  the  vortex  is  increased  until  the 
predicted  instantaneous  distortion  factor,  determined  from  the  graphical 
solution,  is  reached. 

Comparisons  presented  by  Melick,  which  were  based  on  two  sets  of 
inlet  data,  showed  good  agreement  with  measured  results.  No  comparisons 
were  n.ade  between  predicted  and  measured  maximum  pressure  contour  maps. 
The  author  also  examined  the  effect  of  using  fewer  compressor  face  RMS 
turbulence  measurements. 
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SECTION  II 


BACKGROUND 

Several  statistical  methods  have  been  developed  to  predict  the 
maximum  distortion  level,  aitd  in  some  instances,  provide  a  predicted 
pressure  contour  map.  Those  analyses  which  offer  a  predicted  map  ca- 
pabil’ty  are  discussed  here.  Two  methods  that  use  a  form  of  superpo¬ 
sition  are  addressed  first,  followed  by  a  discussion  of  three  methods 
that  utilize  a  random  number  process. 

King,  Schuerman,  and  Muller  (Reference  6)  reported  on  a  distortion 
synthesis  procedure  used  to  estimate  the  maximum  inlet  distortion  level 
associated  with  "drift"  stalls  that  occurred  during  an  inlet-engine  test. 
The  procedure  consists  of  defining  a  reference  180°  segment  (based  on  the 
steady-state  pressure  contour  map)  defined  as  the  segment  with  the 
minimum  steady-state  total  pressure  average.  The  steady-state  pressure 
contour  map  is  intensified  by  subtracting  total  pressure  increments  from 
the  steady-state  probe  pressures  within  the  reference  180°  segment  and  by 
aoding  increments  to  the  probe  pressures  outside  the  reference  segment. 
The  total  pressure  increment  is  defined  as  the  sum  of  the  root-mean-square 
(RMS)  of  the  total  pressure  fluctuations  for  a  given  probe  and  the 
average  RMS  for  all  the  probes  contained  within  the  associated  180' 
segment.  The  pressure  increment  is  multiplied  by  a  coefficient  determinec 
by  the  extreme  value  extrapolation  method  of  Jacocks  and  Kneile  (Reference 
7).  Briefly,  Jacccks'  method  assumes  that  a  short  time  segment  of 
distortion  data  can  be  used  to  statistically  predict  the  maximum  dis¬ 
tortion  level  corresponding  to  any  time  period  of  inlet  operation.  King 
concluded  that  there  v/as  good  agreement  between  the  synthesized  and 
measured  distortion  with  a  tendency  to  overpredict  the  maximum  distortion 
level.  Unfortunately,  no  comparisons  were  made  between  the  intensified 
and  measured  pressure  contour  maps  for  the  maximum  distortion  levels. 

Melick's  analytical  model  (References  8  and  9)  has  a  statistical 
basis  that  is  intended  to  provide  both  an  understanding  and  quantitative 
description  of  unsteady  inlet  flow.  A  solution  to  the  one-dimensional, 
time-dependent,  Navier-Stokes  equation  is  used  to  describe  the  flow  field 
of  an  isolated  vortex.  It  is  hypothesized  that  the  flow  field  consists 


pressure  to  the  steady-state  total  pressure).  In  spite  of  these  results, 
the  approach  of  using  random  numbers  to  generate  a  synthesized  pressure 
distortion  map  remained  promising  if  improvements  could  be  made  to  the 
method. 

In  Section  III,  a  description  of  the  basic  statistical  prediction 
nicdel  is  presented.  The  inlet  Jota  base  and  the  engine  distortion 
parameters  used  to  determine  maximum  distortion  levels  are  described. 
Comparisons  between  the  model  and  DYNADEC  analysis  results  are  presented 
to  show  the  limitations  of  the  method  in  predicting  the  maximum  distor¬ 
tion  level  and  pressure  contour  map  over  a  range  of  inlet  distortion  and 
turbulence  levels.  It  is  generally  assumed  that  inlet  total  pressures 
are  normally  distributed  and  uncorrelated.  Previous  studies  (References 
?.  and  5)  have  investigated  and  supported  the  assumption  of  normality  and 
therefore  will  not  be  addressed.  The  assumption  that  the  pressures  are 
uncorrelated  is  examined  to  assess  its  importance  to  the  method.  Because 
the  power  spectrum  of  the  pressures  used  in  a  distortion  analysis 
exhibits  a  decreasing  amplitude  with  increasing  frequency  due  to  turbu¬ 
lent  mixing  and  filtering  to  define  engine  sensitivity,  an  additional 
constraint  is  imposed  on  the  statistical  model.  Two  digital  filters  are 
described  that  shape  the  power  spectral  density  (PSD)  of  the  random  numbers 
set  to  approximate  the  PSD  of  inlet  pressure  data. 

Section  IV  focuses  on  the  revised  statistical  prediction  model  that 
incorporates  the  two  digital  filters  described  in  Section  III.  Compari¬ 
sons  between  predicted  and  measured  maximum  distortion  levels  and 
pressure  contour  maps  are  presented  to  demonstrate  the  applicability  of 
the  model.  A  map  averaging  approach  is  then  described  that  offers  a 
substantial  improvement  in  the  quality  of  the  pressure  contour  map.  This 
study  also  investigates  the  use  of  a  reduced  number  of  measured  RMS 
turbulence  values  in  the  analysis.  Results  are  presented  based  on  using 
eight  turbulence  values. 

Section  V  presents  a  summary  and  conclusions  of  the  investigation. 


models  in  large  facilities,  and  the  extensive  instrumentation  and  data 
processing  requirements.  Thus  smaller  models,  with  less  instrumentation 
and  less  data  processing  requirements,  are  being  considered  for  use  in 
smaller  facilities.  In  addition,  there  has  been  an  effort  to  develop 
analytical  models  which  are  based  on  a  more  simplified  data  processing 
approach  to  supplement  the  analog  screening  process  described  above. 

A  number  of  technical  needs  can  be  fulfilled  by  having  an  ar.alytical 
model  that  synthesizes  time-variant  distortion.  Such  a  tool  would  be 
useful  in  the  inlet  preliminary  design  and  development  phases.  The 
availability  of  such  a  model  could  provide  on-line  prediction  capability 
during  wind  tunnel  tests,  allow  the  prediction  of  peak  distortion  based 
or  a  limited  amount  of  pressure  instrumeritation ,  end  provide  space  at  the 
compressor  face  plane  to  measure  other  flow  quality  parameters.  Such  a 
method,  however,  must  provide  the  pressure  contour  pattern  as  well  as  the 
maximum  oistortion  level,  be  accurate  at  high-turbulence  levels,  and  be 
economical  relative  to  the  testing  with  full  dynamic  instrumentotion. 

In  recent  years  various  inlet  flow  distortion  analyses  have  beer 
proposed  which  use  statistics.  These  methods  use  either  the  distortion 
parameter  or  the  pressure  statistical  properties  to  predict  the  maximum 
distortion  level.  Some  of  these  methods  include  a  synthesis  of  the 
pressure  distortion  pattern.  Distortion  pattern  synthesis  has  been 
accomplished  by  the  intensification  of  the  steady-state  pattern  and  by 
the  use  of  a  random  number  process.  These  synthesis  methods  are 
described  in  Section  II. 

It  is  the  use  of  the  random  number  process,  coupled  with  the  inlet 
pressure  statistical  properties  that  is  the  subject  of  this  report.  This 
approach  has  been  examined  previously  by  three  investigators  (References 
2,  3,  A).  Because  of  the  potential  of  these  methods,  a  model  similar  to 
that  described  by  Stevens  (Reference  4)  was  used  in  comparison  with 
several  inlet  distortion  data  sets.  The  initial  results  appeared  to 
validate  the  model,  particularly  in  terms  of  the  predicted  peak 
distortion  level.  As  cases  with  increasing  average  inlet  turbulence  were 
investigated  and  as  fewer  RMS  turbulence  mieasuretnents  were  used  in  the 
analysis,  the  quality  of  the  predicted  maps  deteriorated.  (Turbulence  is 
defined  as  the  ratio  of  the  standard  deviation  or  RMS  of  the  tine-variant 
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Figure  7  Basic  Method  (RfO  -  Pressure  Contour  Map  and 
Histogram,  [■loder^te  Turbulence 


-10 


-5 


0 


b 


10 


-1  D 


niFFFPENCE  IN  PROBF  RFmVFRY 

Fiqure  8  Basic  Method  (GP'-,)  -  Pressure  Contour  flap 
Hiqh  Turbulence 


and  Histoqram, 


cases,  the  predicted  maximuni  distortion  (APRSp)  levels  are  greater  than 
the  measured  values  including  the  low-turbulence  case  (T  -  .014),  where 
there  is  good  agreement  between  the  predicted  and  measured  pressure 
contours.  As  turbulence  increases,  there  is  more  disparity  between  the 
predicted  and  measured  distortion  maps.  The  histograms  indicate  a 
decreasing  concentration  of  probes  centered  about  the  zero  difference  in 
predicted  versus  measured  pressure  recovery  with  a  wider  variation  in 
that  difference.  For  example,  only  fifteen  probes  have  predicted  pres¬ 
sures  within  ±2  percent  of  their  measured  values  for  the  high-turbulence 
case  (T  =  .040)  compared  to  twenty-two  probes  within  that  band  for  the 
moderate-turbulence  case  (T  =  .027).  Further,  there  is  a  22  percent 
variation  in  the  pressure  recovery  difference  for  the  high-turbulence 
case  compared  to  a  range  O'f  17  percent  for  the  moderate-turbulence  case. 

The  capability  of  the  basic  model  using  the  PWA  distortion  parame¬ 
ters  is  illustrated  in  Figures  9  through  13.  A  cut-off  filter  frequency 
of  1000  HZ  was  used  to  filter  the  pressure  data.  The  differences  between 
the  predicted  and  measured  total  distortion,  along  with  its  components 

K 

are  shown  in  Figure  9.  The  predicted  values  of  within  a  ±20 

percent  band,  showing  a  tendency  to  overpredict  the  measured  peak  distor¬ 
tion  levels  at  low  distortion  levels  and  underpredict  at  higher  distortion 
levels.  K-  shows  a  similar  trend  since  it  constitutes  a  substantial 
portion  of  the  total  distortion  level.  The  radial  distortion  component, 
'RA^  exhibits  a  substantial  variation  with  differences  as  large  as  50 
percent. 

Figure  10  presents  the  %RMSpjpp  as  a  function  of  average  turbulence. 
As  before,  the  data  points  represent  the  cases  presented  in  the  previous 
figure  while  the  linear  regression  line  is  based  on  the  repetition  of 
several  cases  repeated  several  times  using  different  sets  of  random 
numbers.  The  increased  slope  of  the  regression  line,  compared  to  the 
slope  of  the  regression  line  in  Figure  5,  reflects  a  greater  difference 
between  predicted  and  measured  pressures  and  is  attributed  to  filtering 
the  pressure  data  at  a  higher  cut-off  frequency. 

Comparisons  between  the  predicted  and  DYMAOEC  generated  pressure 
contour  maps  are  shown  in  Figures  11,  12,  and  13.  For  the  low-turbulence 
case.  Figure  11,  there  is  reasonable  agreement  between  the  two  maps  with 
most  cf  the  disparity  in  the  contour  shapes  between  the  eight  and  eleven 
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o'clock  positions.  Thirty-five  of  the  predicted  pressure  recoveries  are 
within  ±2  percent  of  their  measured  values.  The  predicted  contour  maps 
for  the  moderate-  and  high-turbulence  cases.  Figures  12  and  13,  bear 
little  resemblance  to  the  measured  contour  maps,  despite  the  fact  that 
the  predicted  distortion  levels  are  in  good  agreement  with  the  measured 
peak  values.  The  histograms  show  the  same  trend  as  before,  as  turbulence 
increases,  there  is  a  flattening  and  spreading  of  the  distribution  of  the 
difference  between  predicted  and  measured  pressures.  For  the  moderate- 
turbulence  (T  =  .027)  case,  twenty  probes  agree  within  ?2  percent  of 
their  measured  values  with  a  range  of  19  percent  in  the  difference 
between  predicted  and  measured  recoveries.  The  number  of  probes  agreeing 
to  within  ±2  percent  decreased  to  twelve  while  the  range  increased  to  30 
percent  for  the  high-turbulence  case  (T  =  .040). 

Figures  14  through  18  present  a  set  of  comparisons  for  the  second 
set  of  GE  engine  distortion  parameters  (GF^).  Inlet  pressure  data  was 

i 

filtered  at  a  cut-off  filter  frequency  of  500  HZ.  The  differences 
between  predicted  and  measured  values  of  IDL, 

shown  in  Figure  14.  Predicted  values  of  IDL  vary  from  essentially  no 
difference  to  values  in  excess  of  50  percent  from  measured  values.  The 
circumferential  distortion  exhibits  a  similar  behavior  as  it  is  the 
dominate  term  of  the  total  distortion.  The  radial  distortion  is  randomly 
distributed  with  a  wide  variation  from  -25  to  ■'■60  percent  of  the  measured 
values. 

The  average  standard  deviation  of  the  difference  in  probe  recoveries 
for  this  data  set  is  shown  in  Figure  15.  As  before,  the  linear  regres¬ 
sion  line  is  based  on  repeating  several  cases  several  tiriies  to  account 
for  variations  due  to  the  random  number  process.  The  resu1tir,u  slope  of 
the  regression  line  is  similar  to  that  of  the  GF^  distortion  parameter 
set.  Figure  5,  where  the  pressure  data  was  filtered  at  500  HZ. 

A  comparison  of  predicted  and  measured  pressure  contours  are  shown 
in  Figures  16,  17,  and  18.  With  some  exception,  the  predicted  dr,(‘ 
measured  contour  lines  for  the  low-turbulence  case.  Figure  16,  shew  good 
agreement  with  each  other.  Twenty-seven  probes  agree  to  within  ? 
percent  of  their  measured  pressures.  The  range  of  the  difterence  is  ’1 
percent.  For  the  nioderate-turbulence  case,  Figure  17,  the  predicted  and 
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Figure  13  Basic  Method  (PMA)  -  Pressure  Contour  Map  and  Histooran, 
High  Turbulence 
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DIFFERENCE  IN  PROBE  RECOVERY 


i'urf  16  Basic  Method  (GF^)  -  Pressure  Contout 
Histo<jrapi,  low  Turbulence 


of  filter  is 


(10) 


M 

Y,  -  Vn 

k=l 

A  filter  with  three  coefficients  was  selected  to  minimize  model  points  in 
the  power  spectrum.  Generalized  power  spectral  density  shapes  for  a 
random  number  set  filtered  with  the  non-recursive  filter'  are  illustrated 
in  Figure  26.  As  the  values  of  hj^  increase  and  h^  decrease,  the  spectrum 
exhibits  a  progressively  steeper  roll-off  characteristic.  The  influence 
ot  this  filter  on  the  predicted  maximum  distortion  level  and  pressure 
contour  map  is  also  determined  by  the  engine  sensitivity  cut-off  filter 
frequency.  If  the  cut-off  frequency  is  in  the  lower  portion  of  the  spec¬ 
trum,  where  it  is  essentially  flat,  the  effect  of  the  non-recursive 
filter  is  very  small.  Thus,  the  non-recursive  filter  has  a  secondary 
effect  on  the  predicted  maximum  distortion  level  for  small  values  of  h|^ 
and  lower  cut-off  filter  frequencies. 

The  filter  used  to  describe  engine  sensitivity  has  the  greatest 
effect  on  predicted  levels  of  distortion.  Engine  manufacturers  have 
specified  three  or  four-pole  linear-phase  (Bessel)  and  constant-amplitude 
(Butterworth)  filters  to  describe  critical  engine  frequencies  (Reference 
20).  Butterworth  filters  exhibit  a  sharper  roll-off  characteristic 
beyond  the  cut-of'^  frequency  while  Bessel  filters  minimize  in-phase 
relaf ioriships  with  frequency.  Three-pole  Butterworth  filters  are  used  to 
tiller  the  pressure  data  used  in  the  DYNADEC  analysis. 

A  recursive  filter  has  been  developed  to  represent  a  three-pole 
(irrlog  Butterworth  filter.  Whereas  the  output  of  a  non-recursive  filter 
is  a  function  of  the  input  only,  the  recursive  filter  output  is  a  function 
(•^  previous  output  as  well  as  input.  This  recursive  filter  (Reference 
21}  uses  a  bilinear  transformation  of  a  continuous  filter  function, 
defined  as 

S  (11) 

Z+1 

where  S  is  a  complex  variable  and  Z  is  a  rational  function  that  maps  the 
atciginary  axis  of  the  S-plane  onto  the  unit  circle  of  the  Z-plane.  The 
desired  digital  filter  is  obtained  by  substituting  the  bilinear  transform 
into  the  analog  transfer  function,  H(S),  which  for  a  three-pole 
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correlation  is  seen  to  have  an  overall  lower  level  comparea  to  the 
previous  two  figures  due  to  the  relatively  low-turbulence  levels  in  the 
low  pressure  regions. 

These  examples  show  the  complexity  of  describing  the  turbulent  flow 
in  terms  of  a  correlation  coefficient.  Obviously,  it  is  not  appropriate 
to  apply  the  model  to  those  conditions  where  a  discrete  frequency  exists, 
such  as  inlet  buzz,  duct  resonance,  or  any  other  in  or  out-of-phase 
oscillation.  These  conditions  can  not  be  adequately  described  with  only 
the  steady-state  and  RMS  pressures.  A  statistically  significant  level  of 
correlation  (p-->.3)  has  been  shown  to  exist  for  adjacent  probes  on  a 

'  «J 

rake,  in  particular,  those  located  near  the  outer  wall.  In  the  absence 
of  any  flow  oscillations,  it  appears  that  the  assumption  that  the 
pressures  are  uncorrelated  is  not  a  necessary  condition  for  the  model  to 
be  valid. 


Digital  Filter  Development 

The  power  spectral  density  of  inlet  pressure  data  used  in  a  dynamic 
distortion  analysis  has  two  features.  First,  the  spectrum  exhibits  a 
decreasing  amplitude  with  increasing  frequency  as  more  energy  is  usually 
present  at  lower  frequencies.  Second,  the  spectrum  has  a  sharp  roll-off 
above  some  frequency  due  to  the  use  of  a  filter  defining  engine  sensi¬ 
tivity.  Digital  filters  have  therefore  been  added  to  the  basic  prediction 
model  to  impose  the  additional  constraint  that  the  pov/er  spectral  density 
for  the  generated  random  numbers  will  have  a  shape  approximating  that  of 
inlet  pressure  data. 

A  normally  distributed,  uncorrelated  set  of  random  numbers  is 
equivalent  to  white  noise  which  has  a  flat  power  spectrum.  The  PSD 
function  for  low-pass  white  noise  is  defined  as 

(f)  =  a,  0<f<B;  otherwise  zero  (9) 

A  non-recursive  filter  is  used  to  satisfy  the  general  observation 
that  more  energy  is  contained  in  the  lower  portion  of  the  spectrum. 
Non-recursive  filters  represent  a  data  averaging  technique  that  uses 
variable  weighting  coefficients.  The  formula  (Reference  6)  for  this  type 
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cut-off  filter  frequency.  The  more  closely  spaced  probes  near  the  outer 
duct  wall  exhibit  a  high  correlation  level  with  the  level  decreasing 
rapidly  as  probes  become  progressively  further  apart  (probes  located  near 
the  hub  or  two  non-adjacent  probes).  The  correlation  is  also  a  function 
of  cut-off  filter  frequency.  As  cut-off  frequency  decreases,  the 
correlation  coefficient  increases. 

Figure  23  presents  the  correlation  coefficient  for  probes  along 
several  rakes  for  a  circumferential  distortion  pattern.  As  with  the 
results  presented  in  the  previous  figure,  the  level  of  correlation 
decreases  with  increasing  distance  between  probes.  There  is  no  system¬ 
atic  trend  in  the  correlation  level  for  probes  located  progressively 
between  the  high  and  low  pressure  rakes.  In  addition,  the  range  of  the 
correlation  level  for  the  various  rakes  appears  to  be  within  the  scatter 
of  the  correlation  for  the  high-  and  low-pressure  regions  denoted  by  Mace 
and  Sedlock  (Reference  18). 

The  correlation  for  another  circumferential  distortion  pattern  is 
presented  in  Figure  24.  In  addition  to  a  wide  variation  in  the  coeffi¬ 
cient  for  a  given  probe-to-probe  spacing,  the  correlation  exhibits  a 
"sawtooth"  characteristic.  This  "sawtooth"  behavior  illustrates  that  the 
correlation  is  a  function  of  other  variables  as  well  as  distance.  In  a 
study  by  Martin  (Reference  5),  for  example,  the  cross-correlation  and 
phase  angle  showed  significant  differences  between  adjacent  compressor 
face  probes.  To  explain  the  behavior  of  the  cross-correlation  function, 
a  turbulent  zone  concept  was  proposed  to  describe  the  boundary  layer/shock 
wave  mechanisms  present  in  the  flow.  Other  mechanisms  that  would  affect 
the  correlation  function  are  boundary  layer  separation  ana  the  interaction 
between  shock  waves.  For  the  correlation  to  be  only  a  function  of 
distance,  the  turbulent  flow  must  be  isotropic,  or  at  least  locally 
isotropic,  that  is,  the  decay  in  turbulence  has  no  preference  for  any 
specific  direction  (Reference  19).  Such  conditions  can  be  approximated 
behind  a  grid  with  no  pressure  distortion  present.  Thus,  the  "sawtooth" 
behavior  of  the  correlation  is  more  likely  representative  of  inlet 
diffuser  flow. 

Figure  25  presents  the  correlation  for  probes  located  on  rakes  in 
t^he  low  pressure  region  from  three  radial  distortion  patterns.  The 
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correlation  between  probes  located  along  a  rako.  Ihe  cross-correlation 
function  is  obtained  by  averaging  the  product  of  two  random  pressure-time 
histories  over  a  statistically  significant  time  period.  The  cross- 
correlation  function  for  two  laterally  spaced  probes  is 

=  —  /<(t)y(t)dt  (5) 

xy  ^  j 

0 

The  correlation  coefficient  is  defined  as 


it)  =  -  (6) 

APy(x)  AP.p(y) 

RMS  RMS 

r  P 

where  ‘  T^^^^  is  the  standard  deviation  of  a  time-variant  pressure. 
Filtered  time-variant  pressure  data  were  input  to  these  equations  which 
were  programmed  on  ars  analog  computer.  The  data  was  filtered  at  cut-off 
frequencies  of  12S  HZ  X  where  N  =  0,  to  deteiMiine  the  effect 

of  frequency  on  the  RMS  pressure  and  the  correlation  coefficient. 

Figure  21  presents  the  variation  in  RMS  pressure  as  a  function  of 
frequency  for  two  inlet  data  sets.  The  RMS  pressures  for  each  probe  for 
all  the  cases  used  in  the  investigation  are  included  in  this  figure.  The 
RMS  pressures  were  non-dimensional ized  by  the  maximum  RMS  pressure  at  the 
maximum  frequency.  A  logarithmic  transformation  (Reference  15)  of  the 
form 

Y  =  Bo  ^  where  X'  =  logj^^X  (7) 


was  used  to  fit  a  regression  curve  to  the  data.  The  regression  equations 
tor  both  data  sets  are  virtually  identical  and  therefore  described  by  a 
single  expression 


RMS^ 

- -  1.020  +  .520  log^^ 

^^^MAX 


Ferre  ot  the  correlation  results  obtained  are  presentee  in  the 
■ficiurfs  ^hat  follow.  Figure  22  presents  the  correlation  coefficient, 

I  ,  as  a  function  of  the  distance  between  probes  alc'og  a  rake  and 

‘  J 
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Figure  19  Probe- to -'^ro be  Cross  Correlation 


Figure  20  f’robe-to-r’robe  Cross  Correlation 


the  probes  are  more  than  10  percent  of  the  duct  radius  apart.  The 
correlation  on  which  Motycka  based  his  conclusion  is  shown  in  Figure  19 
and  is  based  on  pressure  measurements  taken  behind  a  turbulence  generator 
with  grill,  pipe,  and  egg  crate  devices  positioned  between  the  generator 
and  the  engine  compressor  face  (Reference  16).  The  correlation  was 
developed  to  indicate  the  scale  of  the  turbulence  eddies  generated  by  the 
different  devices.  For  a  typical  spacing  between  probes  (?0  percent  of 
the  duct  radius),  the  correlation  coefficient  has  an  approximate  value  of 
five-tenths.  Motycka  pointed  out  that  a  high  correlation  would  exist 
across  the  engine  face  under  the  conditions  of  inlet  buzz  and  duct 
resonance.  In  those  instances,  the  engine  reaction  to  the  discrete 
frequency  would  be  evaluated  separately. 

Borg  (Reference  2)  also  cautioned  against  applying  the  prediction 
model  to  low  mass  flow  conditions  where  inlet  buzz  may  exist,  but 
presented  no  data  in  support  of  the  assumption  that  the  pressures  are 
uncorrelated. 

Crites  and  Heckart  (Reference  17),  reviewing  model  scale  and 
transducer  spacing  requirements,  discussed  the  need  to  have  the  distance 
between  probes  sufficiently  small  to  allow  an  accurate  interpolation  of 
time-variant  pressure  data.  The  correlation  fields  of  two  probes  would 
therefore  have  to  overlap  at  a  statistically  significant  value,  implied 
to  be  three-tenths. 

Mace  and  Sedlock  (Reference  IS)  explored  the  use  of  a  spatial 
correlation  coefficient  for  estimating  the  time-variant  pressure  at  some 
point  between  two  laterally  spaced  pressure  probes.  The  spatial  corre¬ 
lation,  shown  in  Figure  20,  was  used  in  a  set  of  regression  equations  to 
estimate  the  pressure  fluctuations.  Analysis  showed  that  the  waveform  of 
thf  estimated  pressure  fluctuations  were  representative  of  the  frequency 
and  amplitude  characteristics  of  the  measured  pressure  waveform  at  ihat 
location.  It  was  concluded  that  adjacent  high-response  probes  on  a  rake 
exhibited  a  statistically  significant  cc-rrelation.  comparison  of  the 
correlation  developed  by  Mace  and  Sedlock  with  that  presented  by  Motycka 
show  the  correlations  to  be  of  comparable  magni'*^udp. 

As  part  of  this  investigation,  a  correlation  study  was  acconf)!  i  shed 
with  a  porMon  of  the  inlet  pressure  data  to  determine  the  level  of 
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measured  maps  exhibit  yood  similarity,  but  the  number  of  probes  within  I'd 
percent  of  their  measured  values  has  decreased  to  eighteen  and  the  range 
of  the  difference  has  increased  to  18  percent.  For  the  high-turbulence 
case  (T  =  .047)  shown  in  Figure  18,  there  is  no  similarity  between  the 
predicted  and  measured  contour  maps.  Only  fourteen  probes  are  within  2 


percent  of  their  measured  values  while  the  range  has  increased  to  32 
percent. 

In  summary,  it  appears  that  a  statistical  prediction  model's  ability 
to  accurately  predict  the  maximum  distortion  level  may  have  little  to  do 
with  providing  an  accurate  distortion  map.  A  distortion  parameter  can 
provide  some  averaging  that  minimizes  the  effect  of  the  differences 
between  predicted  ana  measured  probe  recoveries.  For  example,  the  PWA 

circumferential  distortion  parameter,  K.,  uses  a  Fourier  curve  fit  of  the 

0 

pressure  distribution  about  a  ring,  and  consequently  may  be  less  affected 
by  errors  in  pressure  values.  On  the  other  hand,  the  GE  circumferential 
distortion  parameter,  IDC,  uses  the  low  recovery  probe  value  in  defining 
the  distortion  level.  If  that  minimum  value  is  in  error,  it  has  a  sig¬ 
nificant  effect  on  the  prediction  of  total  distortion. 

Histograms  for  the  forty-nine  cases  investigated  have  shown  that  as 

turbulence  increases,  there  is  a  flattening  and  spreading  of  the 

distributior,  with  large  differences  between  predicted  and  Pleasured  probe 
pressures . 

Correlation  Analysis 

An  assumption  of  the  basic  statistical  prediction  model  is  that  the 
individual  compressor  face  pressures  are  uncorrelated.  If  this  assump¬ 
tion  is  true,  then  the  pressure  fluctuations  at  one  location  have  no 
bearing  upon  the  pressure  fluctuations  at  anotfter  location  and  cannot  be 
used  to  predict  anything  about  the  pressures  at  that  point.  However, 
previous  investigatiof'S  (References  5,  17,  and  18)  have  shown  that  a 

statistically  significant  level  of  correlation  can  exist  between  two 

laterally  spaced,  adjacent  probes.  The  assumption  tliat  the  pressures  are 
uncorrelated  was  therefore  examined  to  determine  its  relative  importance 
to  the  model . 

Motycka  (Reference  3)  presented  data  to  shew  that  the  correlation 
between  aojacent  probes  is  generally  low  (less  than  60  percent)  because 
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Figure  I 


Basic  Method  (GEo)  -  Pressure  Contour  Map  and  Histogram, 
Moderate  Turbulence 


Butterworth  filter  is 
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where  is  an  analog  frequency  variable  defined  as 

“‘a  ~  sampling  rate 


(13) 


Fxpanding  the  equation,  and  combining  certain  terms  defined  as: 
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the  following  coefficients  can  then  be  defined: 

^0  ""  ®3 

=  02  “  3(1)^  /Aq 
~  ^2  ~  ^9  ~ 

Then,  the  filter  output  becomes. 


(lA) 


(15) 


Yp  =  BpX.  +  BjZ‘^X.  +  E2Z'^X.  +  B3Z’^X.  +  C^Z’^^  +  C.Z'^  +  C3Z-^Yp(16) 
-1  -2  -3 

where  X^,  Z  X^. ,  Z  X^,  Z  X^  are  the  current  and  past  three  input  values 

-1  -2 

of  X^,  respectively,  and  Z  Y^,  Z  Y^  and  Z  "Y^  are  the  past  three  output 
values. 

By  defining  the  Butterworth  filter  in  this  manner,  the  coefficients 
can  be  determined  by  specifying  the  cut-off  frequency  and  sampling  rate. 
Based  on  experience,  the  sampling  rate  should  be  at  least  four  times  the 
cut-off  filter  frequency  as  lower  sampling  rates  were  found  to  introduce 
instabilities  in  the  filter  algorithm.  Figure  27  presents  an  example  of 
the  smoothed  power  spectral  density  for  a  random  number  set  filtered  at  a 
cut-off  frequency  of  500  HZ. 
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The  problem  is  magnified  in  implementing  these  two  digital  filters 
because  there  are  forty  independent  synthesized  pressure  signals  to  be 
developed.  Rather  than  create  forty  long  strings  of  synthesized  pres¬ 
sures,  a  moving  window  concept  and  several  small  arrays  are  used  to  store 
current  and  past  input  and  output  filter  values.  A  two  step  process  is 
used,  an  initialization  phase  to  start  the  filter  process,  and  a  run 
phase  to  filter  the  synthesized  pressures.  The  following  is  a  brief 
description  of  the  arrays  used  in  the  filter  process: 

Step  Operation 

1  Build  an  array  (RANO)  of  120  random  numbers, 

RANO  =  A(I),  I  =  1,120 

2  Build  an  array  (RAN)  of  forty  random  numbers, 

RAN  =  B(I),  I  =  1,40 

2  Slii-'^t  RANO  one  place  to  the  left  and  replace 

each  third  value  from  the  RAN  array, 

RANO  =  A(2),  A(3),  B(l),  A(5),  A{6), 

B(2),...,  A(119),  A(120),  B(40) 

4  Filter  the  RANO  array  using  the 

non-recursive  filter  and  create  an  array 
(PA40H)  of  forty  values, 

PA40H  =  C(l),  C{2),..,  C{40) 

where:  C(l)  =  hj^  A(2)  +  h^  A(3)  +  h,^  P(l) 

C(2)  =  h^  A(5)  +  h^^  A(6)  +  B(2) 

C(40)  =  h|^  A{119)  +  hp  A(120)  +  h|.  B(40) 

L  Sfn'ft  an  array  (PA120),  initially  filled  with 

zeros,  one  place  to  the  left  and  replace  each 
third  value  with  a  value  from  the  PA40H  array, 

PA120  =  0,  0,  C(l),  0,  0,  C(?) .  0,  0,  C(40) 


Steps  2  through  5  are  repeated  three  times  during  the  initialization 
phase.  The  RAN  array  is  filled  with  a  new  set  of  random  numbers  with 
each  iteration. 


At  the  end  of  the  initialization  phase,  the  PA120  array  has  been 
filled  with  the  output  of  the  PA40H  array  (the  non-recursive  tiltpr 
output).  The  120  numbers  in  the  PA120  array  consist  of  forty  independent 
sets  of  three  numbers  that  represent  three  samples  of  each  synthesized 
pressure.  The  first  two  samples  are  past  non-recursive  filter  values 
while  the  third  sample  is  the  current  value.  The  current  value  is 

determined  from  the  past  two  values  from  RANO  and  a  new  random  number 
from  RAN. 

In  the  run  phase,  steps  2  through  5  are  continued  to  generate  the 
forty  filtered  (non-recursive)  synthesized  pressure  samples.  The 
Butterworth  filters  are  then  activated.  Two  additional  arrays  are 
employed,  an  array  which  holds  the  current  forty  values  and  another  array 
which  contains  forty  sets  of  the  past  three  output  values  of  the 
Butterworth  filters.  By  manipulating  these  arrays  in  the  same  manner  as 
described  above,  the  Butterworth  filter  can  be  programmed  using  a 
relatively  small  amount  of  storage. 
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SECTION  IV 


IMPROVED  MODEL  ANALYSIS 

This  chapter  begins  by  describing  the  improved  statistical  predic¬ 
tion  model  that  includes  the  two  digital  filters  described  in  Section 
III.  Predictions  from  this  model  are  compared  to  the  same  inlet  distor¬ 
tion  data  used  to  examine  the  capabilities  of  the  basic  model.  A  map 
averaging  concept  is  described  which  offers  a  substantial  improvement  in 
the  quality  of  the  predicted  pressure  distortion  maps.  Finally,  the 
results  using  eight  measured  RMS  turbulence  values  in  the  analysis  are 
presented. 


Improved  Statistical  Prediction  Model  Description 

A  schematic  of  the  improved  model  is  shown  in  Figure  28.  The  method 
consists  of  three  elements;  the  generation  of  the  compressor  face  dynamic 
total  pressures,  the  determination  of  the  maximum  distortion  level  and 
pressure  contour  map,  and  the  map  averaging  method  that  provides  the  most 
probable  maximum  distortion  level  and  pressure  contour  map.  A  more 
detailed  discussion  of  map  averaging  is  presented  later. 

The  RMS  turbulence  and  random  numbers  are  combined,  as  was  done  in 
the  basic  method.  However,  these  synthesized  pressure  signals  are  now 
input  into  the  two  digital  filters  added  to  the  basic  model.  The  first 
filter  is  the  non-recursive  filter  that  provides  a  slight  roll-off 
characteristic  over  the  entire  power  spectrum  and  the  second  filter,  the 
recursive  filter,  accounts  for  engine  sensitivity.  The  second  filter's 
output  is  the  filtered  synthesized  fluctuating  pressure  component  that  is 
added  to  the  steady-state  pressure  to  form  the  dynamic  total  pressure. 
The  dynamic  total  pressures  are  input  to  the  distortion  parameter  equation 
where  the  distortion  level  is  determined.  Forty  new  random  numbers  a>-e 
generated,  providing  forty  dynamic  pressures  for  another  equivalent  time 
slice.  The  distortion  level  is  then  computed  and  compared  to  the  current 
maximum  value.  The  larger  value  is  retained,  including  the  pressures  for 
the  distortion  map.  The  sequence  is  repeated  until  the  desired  sample 
size  is  reached  which  is  based  on  cut-off  filter  frequency.  The  solution 
is  restarted  with  another  set  of  random  numbers,  another  maximum 


distortion  map  is  generated,  and  the  sequence  is  repeated  several  times. 
The  pressure  recoveries  of  each  probe  for  all  the  generated  maximum 
distortion  maps  are  averaged  to  develop  the  most  probable  maximum 
distortion  map. 

Prediction  Model  with  Filtering 

The  same  set  of  comparisons  is  used  to  demonstrate  the  capabilities 
of  the  model.  To  reiterate,  the  differences  between  the  predicted  and 
measured  distortion  levels  are  presented  as  a  function  of  measured 
distortion.  Next,  the  standard  deviation  of  the  difference  between 
predicted  and  measured  probe  pressure  recovery  is  examined.  Then,  a 
comparison  is  made  between  the  predicted  and  measured  pressure  contour 
maps,  including  the  histogram  depicting  the  distribution  of  the  differ¬ 
ence  in  probe  recovery.  An  addif-jiial  figure  is  presented  to  further 
illustrate  the  difference  between  the  improved  and  basic  prediction 
models.  A  histogram  is  presented  showing  the  distribution  of  probe 
recovery  difference  based  on  the  average  of  several  solutions  for  the 
particular  example. 

The  improved  model's  capabilities  are  illustrated  first  without  map 
averaging  to  show  the  benefits  of  filtering  the  synthesized  pressures  to 
account  for  engine  sensitivity.  For  the  cases  presented,  the  non¬ 
recursive  filter  is  not  included  in  the  model  for  two  reasons;  first,  the 
cut-off  filter  frequency  (f  =  500  HZ)  used  in  conjunction  with  the  GE 
distortion  parameters  is  in  the  flat  portion  of  the  shaped  spectrum,  and 
second,  the  basic  and  improved  models  with  the  PWA  parameters  tend  to 
underpredict  the  peak  distortion  level  at  higher  distortion  levels. 

Figure  29  presents  the  differences  between  the  predicted  and 
measureo  peak  distortion  levels,  including  the  circumferential  and  radial 
distortion  components,  using  the  GE,  engine  distortion  parameters.  The 
synthesized  pressures  have  been  filtered  at  a  cut-off  filter  frequency  of 
500  HZ.  Recalling  that  the  basic  method  substantially  overpredicted  the 
peak  APRSp  (Figure  4),  including  the  circumferential  and  radial  compo¬ 
nents,  the  results  from  the  improved  model  with  filtering  show 
significantly  better  agreement  with  all  the  measured  components  of  the 
Distortion.  The  average  difference  for  APRSp  is  approximately  -2.1 


percent  with  a  standard  deviation  of  about  7.5  percent.  The  average 
difference  and  standard  deviation  for  and  -3.5  and  4.3 

percent,  and  1.4  and  11.6  percent,  respectively. 

The  improvement  to  the  model  is  further  emphasized  in  the  difference 
between  the  regression  lines  for  the  basic  and  improved  models  shown  ir’ 
Figure  30.  The  data  points,  shown  in  the  figure,  are  the  same  cases 
presented  in  the  previous  figure  with  the  regression  lines  also  based  on 
repeating  several  of  those  cases  several  times  to  account  for  variations 
due  to  the  random  number  process.  Filtering  the  synthesized  pressures 
has  reduced  the  percent  RMS  difference  in  probe  pressure  recovery  from  35 
to  40  percent  at  the  high-  and  low-turbulence  levels,  respectively. 

Figures  31,  33,  and  33  present  the  predicted  and  measured  pressure 
contour  maps  for  the  low,  moderate,  and  high-turbulence  cases  (GE^) 
presented  earlier.  For  the  low-turbulence  case.  Figure  31,  there  is  an 
improved  definition  of  the  predicted  pressure  contours.  There  is  excel¬ 
lent  agreement  betweer;  the  predicted  and  measured  distortion  compcnents. 
The  percent  RMS  difference  has  been  reduced  by  30  percent  to  a  value  of 
1.20.  A  more  dramatic  improvement  between  the  predicted  and  measured 
distortion  maps  for  the  moderate-turbulence  case  can  be  seen  in  Figure 
32.  Substantial  agreement  exists  between  the  pressure  contours  and 
between  the  predicted  and  measured  distortion  levels.  The  histogram 
shows  a  tighter  distribution  about  the  zero  difference  in  probe  recovery 
with  the  range  of  the  difference  having  decreased  from  17  (basic)  to  12 
percent  (improved).  The  RMS  difference  is  ?.63%,  a  reduction  amounting 
to  21  percent.  Similar  results  ere  achieved  for  the  high-turbulence 

case,  Figure  33.  With  the  exception  of  the  recovery  contours  of  .95,  the 
predicted  and  measured  contours  show  very  good  agreemtoit.  Predicted  and 
measured  distortion  levels  are  identical.  The  7RMSpjpp  decreased  from 
4.60  percent  to  3.34  percent,  the  range  decreased  from  22  to  16  percent, 
and  the  number  of  predicted  pressures  agreeing  within  22  percent  of  their 
nrasured  values  increased  from  fifteen  to  twenty-two  in  number,  compared 
to  basic  model  results. 

Figure  34  presents  the  average  distributions  of  probe  recovery 
difference  of  several  solutions  for  the  three  cases  just  presented.  The 
results  show  a  substantial  iriiprovement  to  the  model  with  filtering. 
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Figure  31  Improved  Method  (GF.)  -  Pressure  Contour  Map 
and  Histogram,  Low  Turbulence  (Filtered) 
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Figure  32  Improved  Method  (GF,)  -  Pressure  contour  Map 
and  Histogram,  Moderate  Turbulence  (Filtered) 
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Figure  4b  Improved  Method  (GE2)  -  Pressure  Contour  Map 
and  Histogram,  High^Turbul ence  (Filtered) 
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Fiqure  43  Improved  Method  (^Fp)  -  •Pressure  fontour 

Map  and  llistoqran,  lov/  Turbulence  (Filtered) 


average  differences  amounting  to  -1.9  percent  and  a  standard  deviation  of 
24.7  percent.  While  the  standard  deviation  appears  large,  the  levels  of 
radial  distortion  are  small  emphasizing  the  relative  differences  between 
the  predicted  and  measured  values  of  IDR^^^. 

A  comparison  of  the  percent  RMS  difference  in  probe  pressure 
recovery  for  the  improved  and  basic  methods  is  shown  in  Figure  42.  As 
with  the  other  distortion  parameters,  filtering  has  a  very  beneficial 
effect  of  reducing  the  difference  between  predicted  and  measured 
pressures.  The  reduction  in  between  the  two  models  varies  from 

no  difference  at  very  low-turbulence  levels  to  approximately  25  percent 
at  high-turbulence  levels. 

Pressure  contour  maps  and  histograms  for  the  three  levels  of  turbu¬ 
lence  previously  presented  are  shown  in  Figures  43,  44,  and  45.  For  the 
low-turbulence  case  presented  in  Figure  43,  there  is  excellent  agreement 
t'etween  the  predicted  and  measured  distortion  maps.  The  predicted  level 
of  IDL  is  approximately  10  percent  greater  than  the  measured  peak  IDL, 
however,  the  differences  between  the  predicted  and  measured  circumferen¬ 
tial  and  radial  distortion  components  are  very  small.  It  is  how  those 
elements  are  combined  that  contributes  to  the  difference  between  the 
predicted  and  measured  value  of  IDL.  Thirty-three  predicted  pressures 
are  within  ±2  percent  of  their  measured  values  and  the  ^RMSpjPP  is  1.72, 
a  decrease  of  12  percent  over  the  results  with  the  basic  method.  The 
improvement  to  the  predicted  distortion  map  for  the  moderate-turbulence 
case,  Figure  44,  appears  to  be  only  somewhat  improved  over  the  results 
attained  with  the  basic  model.  However,  an  examination  of  the  histogram 
shows  that  the  range  of  the  difference  in  probe  pressure  recovery 
decreased  from  22  to  15  percent  with  a  resultant  decrease  in  %RMSpjpp  of 
approximately  32  percent.  The  high-turbulence  case  is  shown  in  Figure 
45.  Again,  there  is  an  improvement  over  what  was  predicted  with  the 
basic  model,  attaining  substantial  improvements  in  the  predicted  dis- 
tortiori  levels,  the  reduction  in  the  range  of  the  difference  and  percent 
PMS  difference  in  probe  recovery.  However,  as  was  the  case  with  the  PWA 
high-turbulence  example,  there  appears  to  be  an  upper  limit  for  applying 
the  improved  model  with  filtering. 

The  average  distribution  of  the  difference  in  probe  recovery  is 
presented  in  Figure  46.  The  distributions  for  the  three  turbulence 
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Fiqure  40  nistribution  of  Probe  Recovery  Difference  (PWA) 


Figure  3d  Improved  Method  (PWA)  -  Pressure  rontour  Map 
and  Histogram,  Moderate  Turbulence  (Filtered) 


distortion  map  for  the  fiioderate-turbulence  case.  Figure  38,  shows  very 
good  agreement  with  the  measured  distortion  map.  Other  than  the  obvious 
difference  in  the  nine  to  eleven  o'clock  position,  there  is  good  agree- 
tiient  between  predicted  and  measured  contours.  While  the  number  of  probe 
recoveries  agreeing  within  ±2  percent  is  Just  slightly  greater  than  that 
for  the  results  with  the  basic  mc-del ,  there  is  a  substantial  reduction  in 
%RMSgjPP  (30  percent)  reflecting  a  narrower  range  in  the  difference  in 
probe  recovery  (U  versus  17  percent).  Finally  the  high-turbulence  case 
(T  =  .040}  is  presented  in  Figure  39.  /ilthough  the  predicted  and  measured 
distortion  maps  are  in  better  agreement,  compared  to  the  basic  model 
results,  there  are  nevertheless  significant  differences  between  the  two 
m.aps.  A  reduction  in  of  28  percent  was  attained  over  the  basic 

model,  but  the  reiudvely  small  number  of  predicted  pressures  agreeing 
within  ±2  percent  of  their  measured  values  obviously  is  a  factor  in  the 
disparity  between  the  predicted  and  measured  maps.  Thus,  the  combination 
of  high  average  value  of  turbulence  (T  =  .040)  and  high  cut-off  frequency 
(1000  HZ)  may  represent  an  upper  limit  in  applying  the  improved  model. 

The  average  distributions  of  the  difference  in  probe  pressure 
recovery  for  these  three  cases  are  shown  in  Figure  40.  The  most  notable 
gains  a  e  at  the  low-  and  moderate-turbulence  levels  with  only  small 
improvements  in  the  distribution  for  the  high-turbulence  case.  A  compar¬ 
ison  oT  this  distribution  with  the  distribution  for  the  fiigh-turbiilence 
case  presented  in  Figure  34  shows  the  relatively  small  number  of  pressures 
within  2  percent  o-*"  their  measured  values. 

Figures  41  through  46  illustrate  the  improved  method's  capabilities 
Witti  the  GE^  distortion  parameters.  The  synthesized  pressures  were 
tillered  at  a  cut-ofi  tilter  frequency  500  HZ.  The  differences 
tctween  predicted  and  measured  peak  distortion  levels  are  shown  in  Figure 
ai.  In  contrast  to  the  basic  method  which  predicted  total  and  circumfer¬ 
ential  aistortion  levels  significantly  greater  than  the  measured  values, 
the  improved  method  predicts  levels  that  are  in  substantial  agreement 
with  the  r,6'()r.urpd  values.  The  average  difference  between  predicted  and 
measured  values  of  IDE  is  -5,1  percent  with  a  standard  deviation  of  11.3 
percent.  Fcr  the  average  difference  is  -.1  percent  with  a 

standard  deviation  of  16.5  percent.  The  large  variation  in  the  radial 
distortion  component  has  been  substantially  reduced  as  well  with  the 


Filtering  the  synthesized  pressure  reduces  the  roncje  cf  the  difference 
between  predicted  and  measured  probe  recoveries  and  concentrates  more 
values  about  the  zero  difference  in  recovery. 

The  results  for  tiie  prediction  model  wuti  the  PWA  distortion  parame¬ 
ters  are  summarized  in  Figures  35  through  40.  The  synthesized  pressures 
were  filtered  at  1000  HZ.  Differences  between  predicted  and  measured 
peak  distortion  levels  are  shown  in  Figure  35.  The  overall  trend  for  the 
three  distortion  parameters  is  similar,  although  somewhat  lower  in  level, 
with  that  obtained  with  the  basic  method  (Figure  9),  The  average  differ¬ 
ence  between  predicted  ana  measured  values  of  is  -8.9  percent  with  a 

L 

standard  deviation  of  11.4  percent.  For  the  average  difference  is 

-7.3  percent  and  the  standard  deviation  of  14.5  percent.  The  average 

difference  and  standard  deviation  between  predicted  and  measured  values 
|( 

of  RA.,  are  -13.5  and  35.7  percent,  respectively. 

L 

The  similarity  in  the  results  suggests  that  the  filtering  of  the 
synthesized  pressures  has  a  much  smaller  influence  on  the  predicted 
pressures  as  compared  to  results  obtained  with  the  GEj  parameters  where 
the  pressures  were  filtered  at  5CC  HZ,  The  modest  reduction  in  the 
percent  RMS  difference  in  probe  recovery,  reflected  in  the  small  differ¬ 
ence  between  the  two  regression  lines  in  Figure  36,  would  tend  to  support 
this  observation.  However,  the  regression  lines  are  based  on  a  greater 
population  (several  cases  were  repeated  several  times)  than  shown  here 
and  the  slope  of  the  line  for  the  filtered  results  is  effected  by  large 
^RMSpypp  values  for  turbulence  levels  exceeding  .06.  If  that  regressior 
line  was  based  solely  on  those  cases  having  average  turbulence  levels  of 
less  then  .06,  the  slope  of  the  line  would  be  lower  ana  a  greater 
improvement  noted. 

Evidence  of  this  o(r:clusion  car.  be  seen  with  the  substantial 
improvement  T(>  the  pressure  contour  naps  and  the  distribution  of  the 
difference  in  probe  pressure  recovery  presented  in  Figures  37,  38,  and 
39.  For  the  low-turbulence  case.  Figure  37,  there  is  excellent  agreement 
between  pressure  contours  of  t  [le  predicted  and  measured  distortion  maps 
with  thirty-five  predicted  pressures  within  3  percent  of  their  measured 
values.  The  "PMSj^ypp  is  1.18  which  represents  a  3F  percent  improvement 
(iver  the  result  obtained  with  the  basic  model.  The  predicted  pressure 


levels  show  similar  results  as  noted  previously,  that  is,  most  improve¬ 
ment  occurs  at  low  and  moderate-turbulence  levels.  As  turbulence 
increases,  the  distributions  characteristically  become  flatter  and 
broader,  with  less  than  one-third  of  the  predicted  pressures  being  within 
±2  percent  of  their  measured  values  for  the  high-turbulence  case. 

In  suimary,  it  has  been  shown  that  the  filtering  of  the  synthesized 
pressures  provides  a  significant  improvement  in  the  predicted  distortion 
levels  end  pressure  distortion  maps.  The  quality  of  the  predicted  pres¬ 
sure  distortion  maps  is  a  function  of  the  turbulence  level  and  cut-off 
filter  frequency.  As  the  levels  of  these  two  parameters  increase, 
the  quality  of  the  distortion  map  decreases.  The  upper  limits  for 
applying  the  improved  model  with  filtering  appears  to  be  approximately 
.04  for  turbulence  and  a  cut-off  filter  frequency  slightly  greater  than 
lOOC  HZ. 

Prediction  Model  with  Filtering  and  Map  Averaging 

A  map  averaging  scheme  has  been  included  in  the  model  that  provides  a 
significant  improvement  in  the  accuracy  of  the  pressure  distortion  map. 
The  predicted  distortion  map  and  the  level  of  the  distortion  parameters 
are  the  result  of  a  particular  set  of  random  numbers.  Another  set  of 
random  numbers,  via  a  different  starting  seed  value,  will  generate 
somewhat  different  dynamic  total  pressures  and  hence,  a  different  distor¬ 
tion  nap  and  distortion  parameter  values.  That  difference  is  small  for 
probes  with  low-turbulence  levels,  but  can  be  significant  for  probes 
located  in  regions  of  high-turbulence.  All  the  possible  predicted 
distortion  maps  are  equally  valid,  and  thus  each  predicted  probe  pressure 
is  considered  to  be  part  of  a  distribution.  The  average  maximum  distor¬ 
tion  map,  which  will  be  Defined  as  the  most  probable  maximum  distortion 
map,  is  determined  by  repeating  the  solution  several  times  with  different 
sets  of  random  numbers.  Individual  probe  pressures  are  summed  from  the 
iridividual  distortion  maps  and  averaged  so  that  the  resulting  forty 
pressures  represent  the  data  for  the  most  probable  maximum  pressure 
distortion  map.  For  the  cases  presented,  the  pressures  for  six  maximuni 
distortion  maps  were  summed  and  averaged. 

The  three  previous  cases  with  the  GE^  distortion  parameters  are 
presented  in  Figures  47,  48,  and  49.  Pressure  contours  for  the  predicted 
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map  are  virtually  identical  with  the  measured  map  for  the  low-turbulence 
case.  Figure  47.  Excellent  agreement  exists  between  predicted  and 
measured  distortion  levels.  The  %RMSpjPP  is  .68,  a  reduction  of  60 
percent  compared  to  the  results  attained  with  the  basic  model.  Thirty- 
nine  predicted  pressures  are  within  ±2  percent  of  their  measured  values. 
The  moderate-turbulence  case  is  presented  in  Figure  48.  Again  there  is 
excellent  agreement  between  predicted  and  measured  pressure  contour 
lines.  There  is  correspondingly,  excellent  agreement  between  predicted 
and  measured  distortion  parameter  values.  Twenty-nine  probe  pressures 
agree  to  within  ±2  percent  of  their  measured  values.  The  percent  RMS 
difference  has  decreased  to  1.87,  representing  a  44  percent  improvement 
over  the  basic  model.  The  capabilities  of  the  improved  model  with  the 
map  averaging  are  further  demonstrated  with  the  high-turbulence  case 
presented  in  Figure  49.  The  pressure  contours  of  the  predicted  map  show 
excellent  agreement  with  the  measured  map  with  the  exception  of  the 
pressure  recovery  contour  of  .95.  Predicted  and  measured  distortion 
levels  exhibit  excellent  agreement.  The  number  of  probes  within  ±2 
percent  of  their  measured  values  is  twenty-six,  a  73  percent  improvement 
over  the  results  with  the  basic  method.  The  range  of  the  difference 
decreased  56  percent  and  the  decreased  46  percent  compared  to 

the  basic  model . 

Because  of  the  similar  substantial  agreement  between  predicted  and 
measured  distortion  maps  for  the  low-turbulence  examples,  only  moderate- 
and  high-turbulence  cases  are  presented  for  the  PWA  and  GE^  distortion 
parameters.  Figure  50  presents  the  moderate-turbulence  case  with  the  PWA 
distortion  parameters.  Pressure  contours  for  the  predicted  and  measured 
maps  are  almost  identical  with  the  only  apparent  difference  at  approxi¬ 
mately  the  eleven  o'clock  position.  Predicted  distortion  levels  are  in 
excellent  agreement  with  measured  values.  The  distribution  of  probe 
pressure  recovery  difference  shows  that  twenty-nine  probes  are  within  ±2 
percent  of  their  measured  values.  The  range  of  difference  decreased 
slightly  to  11  percent  with  the  percent  RMS  difference  reduced  to  2.14,  a 
reduction  of  44  percent  over  the  basic  model.  The  high-turbulence  case 
is  shown  in  Figure  51,  For  this  case,  there  is  no  improvement  in  the 
predicted  distortion  map  compared  to  what  was  predicted  with  filtering. 
While  the  range  in  the  probe  recovery  difference  has  decreased,  the 
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Figure  47  Improved  Method  (GEi)  -  Pressure  contour  Map 
and  Histogram,  Low  Turbulence  (Filtered  and 
Map  Averaged) 
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%  DIFFERENCE  IN  PROBE  RECOVERY 
Figure  50  Improved  Method  (PWA)  -  Pressure  Contour 


Map  and  Histogram,  Moderate  Turbulence 
(Filtered  and  flap  Averaged) 


percent  RMS  difference  increased  slightly  over  the  improved  model  with 
filtering.  This  case  has  a  rather  complex  pattern  with  two  high  pressure 
lobes,  and  although  better  results  may  be  achieved  with  a  180°  circumfer¬ 
ential  distortion  pattern,  the  results  again  suggest  an  upper  limit  of 
applicability  of  the  improved  model  with  respect  to  the  combination  of 
cut-off  filter  frequency  (1000  HZ)  and  high-turbulence  level  (.040).  To 
further  define  that  upper  bound,  a  different  case  with  a  turbulence  level 
between  the  moderate-  and  high-turbulence  level  case  is  presented  in 
Figure  52.  Very  good  agreement  exists  between  predicted  and  measured 
contou'^s  with  the  exception  in  the  ten  o'clock  region  of  the  map.  The 
distortion  parameters  are  in  excellent  agreement.  Twenty-four  predicted 
pressures  are  within  ±2  percent  of  their  measured  values. 

The  GE2  moderate-turbulence  level  case  is  presented  in  Figure  53. 
Very  good  agreement  exists  between  predicted  and  measured  pressure 
contours  wit!^  the  exception  of  contours  in  the  four  o'clock  region. 
Excellent  agreement  exists  between  predicted  and  measured  levels  of  total 
and  circumferential  distortion.  The  range  of  the  difference  and  the 
percent  RMS  difference  decreased  28  and  40  percent,  respectively,  over 
the  basic  model.  The  high-turbulence  case  is  shown  in  Figure  54. 
Although  there  is  a  substantial  improvement  in  the  predicted  versus 
measured  contours,  the  distortion  parameter  levels,  the  range  of  differ¬ 
ence  and  percent  RMS  difference  compared  to  the  results  obtained  with  the 
basic  model,  there  nevertheless  are  significant  differences  between  the 
predicted  and  measured  distortion  maps.  Again,  the  applicability  of  the 
improved  model,  including  map  averaging,  in  this  combined  region  of  high- 
turbulence  (T  =  .047)  and  cut-off  frequency  is  questionable.  A  lower 
turbulence  level  case  (T  =  .039)  is  presented  in  Figure  55.  Better 
agreement  exists  between  predicted  and  measured  pressure  contours, 
although  not  as  good  as  that  attained  for  the  GE,  high-turbulence  case. 
Nevertheless,  there  is  excellent  agreement  between  predicted  and  measured 
levels  of  distortion.  Eighteen  probe  pressures  are  within  ±2  percent  of 
their  measured  values. 

The  improvement  to  the  model  with  filtering  and  map  averaging  is 
summarized  in  the  next  six  figures  using  two  measures  of  goodness,  the 
percent  RMS  difference  between  predicted  and  measured  probe  pressure 
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Figure  5?  Improved  Method  (PWA)  -  Pressure  Contour 
Map  and  Histogram,  Moderate  Turbulence 
(Filtered  and  Map  Averaging) 
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Figure  54  Improved  Method  (GE2)  -  Pressure  Contour  Map  and 

Histogram,  High  Turbulence  (Filtered  and  Map  Averaged) 
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Figure  55  Improved  Method  (G'^.p)  -  Pressure  Contour 

Map  and  Histogram,  High  Turbulence  (Filtering 
and  Map  Averaged) 


recoveries,  and  the  number  of  probes  with  predicted  values  within  ±2 
percent  of  their  measured  values. 

The  linear  regression  lines,  Figure  56,  illustrate  the  reduction  in 
the  average  ^RMS^^PP  using  filtering  and  map  averaging  with  the  6E^ 
distortion  parameter  methodology.  With  the  exception  at  low-turbulence 
levels  where  there  is  little  difference  between  the  models,  the  filtering 
of  the  random  numbers,  particularly  in  combination  with  map  averaging, 
provides  a  substantial  reduction  in  the  percent  RMS  difference  between 
predicted  and  measured  pressures.  For  example,  at  an  average  compressor 
face  turbulence  level  of  .04,  the  average  %RMS  difference  with  map 
averaging  and  filtering,  is  2.5,  a  reduction  of  56  percent  compared  to 
the  basic  model. 

The  regression  curves  in  Figure  57  present  the  other  indicator  of 
goodness,  the  number  of  synthesized  pressures  agreeing  within  +2  percent 
of  their  measured  values  using  the  GE^  methodology.  There  is  little 
diftprence  between  the  basic  and  improved  models  for  turbulence  levels  of 
less  than  .01  and  for  levels  greater  than  .06.  Maximum  benefits  with  the 
improved  model  are  in  terms  of  N(±2%)  are  in  .02  to  .04  turbulence  range. 
For  a  turbulence  level  of  .03,  the  combined  filtering  and  map  averaging 
otters  a  substantial  improvement  in  the  number  of  probes  predicted  to  be 
within  -“r  percent  of  their  measured  values.  An  average  of  thirty-two 
pressures  are  within  ±2  percent  of  their  measured  values  for  the  improved 
model,  compared  to  only  eighteen  pressures  for  the  basic  method. 

Figures  58  and  59  present  the  same  summaries  with  the  Pratt 
Whitney  distortion  parameters.  The  reauction  in  average  "RMSpjpp,  Figure 
58,  is  not  as  dramatic  as  it  was  witli  the  GF^  parameters  and  is  primarily 
the  consequence  of  filtering  the  data  at  a  higher  cut-off  filter  fre¬ 
quency.  Nevertheless,  the  filtering  and  map  averaging  still  provide 
significant  improvements  to  the  pressure  contour  maps.  As  was  the  case 
with  the  GE^  parameters,  there  is  little  difference  between  the  basic  and 
improved  models  at  low-turbulence  levels.  For  a  turbulence  level  of  .04, 
the  recuction  in  '^RMSgjPP  amounts  to  33  percent  compared  to  the  basic 
model . 

The  effect  of  the  higher  cut-off  filter  frequency  is  also  seen  in 
Figure  59  for  the  number  of  predicted  pressures  within  i2  percent  of 
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their  measured  values.  The  overall  levels  of  the  regression  curves  with 
the  filtering  and  map  averaging  are  lower  compared  to  the  results  in 
Figure  57.  For  example,  at  a  turbulence  level  of  .03,  an  average  of 
twenty-two  probes  are  within  ±2  percent  of  their  measured  values.  As 
before,  there  is  little  difference  between  the  models  at  higher 
turbulence  levels  (T>.04).  Substantial  benefits  can  be  had  with  the 
improved  method  for  turbulence  levels  as  low  as  .01. 

Summaries  for  the  GE2  distortion  parameters  are  presented  in  Figures 
60  and  61.  The  average  percent  RMS  difference  as  a  function  of 

turbulence  is  presented  in  Figure  60.  As  was  the  case  with  the  other  two 
sets  of  distortion  parameters,  there  is  little  difference  between  the 
basic  and  improved  models  for  turbulence  levels  of  less  than  .01.  The 
reduction  in  the  average  XRMSg^pp  for  the  improved  model  with  map 
averaging,  as  compared  to  the  basic  model  results,  is  not  great  is  that 
attained  with  the  GE^  distortion  parameters.  At  a  turbulence  level  of 
.04,  for  example,  the  reduction  in  %RMSpjpp  amounts  to  39  percent  as 
compared  to  56  percent  with  the  GEj  parameters. 

The  number  of  predicted  pressures  agreeing  within  +2  percent  of 
their  measured  values  with  the  GE^  parameters  is  shown  in  Figure  61.  The 
shapes  of  the  curves  are  the  same  as  those  for  the  GE^  parameters,  but  at 
a  somewhat  lower  level.  At  a  turbulence  level  of  .03,  the  average  number 
of  predicted  pressures  within  ±2  percent  of  their  measured  values  is 
twenty-four.  There  is  little  difference  between  the  models  at  turbulence 
levels  of  less  than  .01  and  for  levels  exceeding  .05. 

In  summary,  it  has  been  shown  that  substantial  improvements  can  be 
made  to  the  predicted  pressure  distortion  maps  by  filtering  the  synthe¬ 
sized  pressure  data  and  by  averaging  the  pressures  from  several  maximum 
distortion  maps.  The  quality  of  the  distortion  map  is  a  function  of  the 
turbulence  level  and  cut-off  filter  frequency  defining  engine  sensi¬ 
tivity.  Relative  differences  between  Figures  56,  58,  and  60  suggest  that 
depending  on  the  inlet  turbulence  level,  the  improved  model  appears 
limited  to  cut-off  filter  frequencies  of  less  than  1500  HZ. 

Prediction  Model  Using  Eight  Turbulence  Values 

The  analysis  thus  far  has  been  based  on  forty  measured  turbulence 
values.  Comparisons  between  predicted  and  measured  pressure  contour  maps 
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are  now  presented  using  eight  measured  values.  The  turbulence 
distribution  for  each  ring  is  assumed  to  be  the  same  as  that  for  the 
third  ring  of  the  eight-rake  by  five-ring  array. 

Contour  maps  and  histograms  based  on  the  GEj  distortion  parameter 
methodology  are  presented  in  the  next  three  figures.  The  low-turbulence 
case,  previously  shown,  is  presented  in  Figure  62.  Compared  to  the 
results  based  on  forty  measurements.  Figure  47,  the  range  of  di^f^ference 
increased  from  4  to  8  percent,  the  number  of  predicted  pressures  within 
±2  percent  of  their  measured  values  decreased  from  thirty-nine  to  thirty- 
five,  and  the  %RMSpjpp  almost  doubled  to  a  value  of  1.30.  Nevertheless, 
there  is  excellent  agreement  between  the  predicted  and  measured  pressure 
contours.  The  mooerate-turbulence  case  is  shown  in  Figure  63.  The  only 
difference  between  the  results  presented  here  and  those  based  on  forty 
measurements.  Figure  48,  is  the  slightly  different  distribution  in  the 
histogram  depicting  the  difference  between  predicted  and  measured  probe 
recoveries.  Otherwise,  there  is  excellent  agreement  between  predicted 
and  measured  pressure  contours  and  the  distortion  parameter  components, 
with  no  change  in  the  range  and  the  number  of  pressures  within  ±2  percent 
of  their  measured  values.  Essentially  the  same  results  were  obtained  for 
the  high-turbulence  case,  Figure  64.  All  the  indications  of  goodness; 
the  pressure  contours,  the  distortion  levels,  and  the  distribution  of  the 
difference  in  probe  recovery  show  virtually  identical  results  to  that 
based  on  forty  turbulence  measurements. 

Figure  65  presents  for  the  moderate-turbulence  case,  the  predicted 
and  measured  contour  maps  with  the  PWA  distortion  parameter  methodology. 
The  agreement  between  the  two  maps  is  considered  very  good,  with  some 
discrepancies  in  the  high  recovery  contours.  The  predicted  and  measured 
distortion  levels  are  in  excellent  agreement.  In  addition,  the  percent 
RMS  difference  between  predicted  and  measured  pressure  recoveries  is 
identical  with  the  results  based  on  forty  measurements.  The  number  of 
probes  agreeing  to  within  ±2  percent  of  their  measured  values  decreased 
to  twenty-six,  compared  to  twenty-nine  based  on  forty  turbulence 
measurements. 

Findly,  the  moderate-turbulence  case  with  the  GEp  distortion 
parameters  is  shown  in  Figure  66.  Virtually  the  same  results  have  been 
attained  as  that  for  the  distortion  map  based  on  forty  measurements, 
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Figure  62  Improved  Method  (GE^)  -  Pressure  Contour 
Map  and  Histogram  -  8  RMS  Turbulence 
Measurements.  Low  Turbulence 
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Figure  63  Improved  Method  (GE, )  -  Pressure  Contour 
Map  and  Histogram  -  8  RMS  Turbulence 
Measurements,  Moderate  Turbulence 
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”1  DIFFERENCE  IN  PROBE  RECOVERY 

Figure  64  Improved  Method  (GE,)  -  Pressure  Contour 
Map  and  Histogram  -  8  RMS  Turbulence 
Measurements,  Hinh  Turbulence 
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Figure  66  Improved  Method  (GEp)  -  Pressure  Contour 
Map  and  Histogram  -‘^8  RMS  Turbulence 
Measurements,  Moderate  Turbulence 
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Figure  53.  The  percent  RMS  difference  is  3.2,  an  increase  of  15  percent. 
The  range  of  difference  has  increased  slightly  and  the  number  of 
predicted  pressures  within  ±2  percent  of  their  measured  value  remains 
essentially  the  same. 

Figure  67  presents  a  comparison  of  %RMSpjpp  based  on  eight  and  forty 
turbulence  measurements  for  all  three  sets  of  engine  distortion  parame¬ 
ters.  With  the  exception  of  the  small  difference  at  low-turbulence 
levels  with  the  and  GE2  distortion  parameters,  the  effect  of  using 
fewer  measured  turbulence  values  on  the  average  percent  RMS  difference 
between  predicted  and  measured  probe  recoveries  is  small.  A  similar 
conclusion  can  be  reached  upon  examining  the  effect  of  fewer  turbulence 
measurements  on  the  average  number  of  predicted  pressures  within  ±2 
percent  of  their  measured  values  shown  in  Figure  68. 

An  examination  of  the  magnitude  of  the  terms  of  the  equation 
defining  dynamic  total  pressure  recovery,  in  particular  the  fluctuating 
pressure  component  term,  would  suggest  that  the  recovery  term  and  the 
random  number  dominate  the  turbulence  value.  Consequently,  only  a 
reasonable  estimate  of  probe  turbulence  is  required  to  predict  the  most 
probable  maximum  distortion  pressure  contour  map.  Previous  investigators 
have  shown  that  the  maximum  distortion  level  can  be  predicted  just  as 
well  with  fewer  turbulence  measurements  as  with  forty  values.  This  study 
has  shown  that  a  good  estimate  of  the  pressure  contour  map  is  possible 
with  fewer  turbulence  measurements  as  well. 
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Figure  68  Predicted  Pressures  within  +2%  Measured  P 


SECTION  V 


SUMMARY  AND  CONCLUSIONS 

A  prediction  method  that  uses  the  inlet  total  pressure  statistical 
properties  and  a  randoni  number  generator  has  been  developed  to  predict 
the  most  probable  maximum  distortion  and  pressure  contour  map.  The  model 
incorporates  two  digital  filters  to  shape  the  random  number  power  spectrum 
to  approximate  that  of  inlet  pressure  data  and  a  map  averaging  approach 
that  substantially  reduces  the  difference  between  predicted  and  measured 
pressures. 

Several  measures  of  goodness  have  been  used  to  assess  the  capability 
of  the  model.  Those  factors  include:  the  predicted  pressure  contour  map 
and  distortion  level,  the  distribution  of  the  difference  and  the  standard 
deviation  of  the  difference  between  predicted  and  measured  pressures,  and 
the  number  of  pressures  agreeing  to  within  +2  percent  of  their  measured 
values. 

The  improved  model  with  filtering  and  map  averaging  offers  a 
substantial  improvement  over  the  basic  model.  Depending  on  the  particu¬ 
lar  engine  distortion  parameter  methodology  and  the  corresponding  cut-off 
filter  frequency  defining  engine  sensitivity,  reductions  in  the  percent 
RMS  difference  between  predicted  and  measured  probe  pressures  varied  from 
30  to  60  percent  at  an  average  inlet  turbulence  level  of  .03.  Corre¬ 
spondingly,  the  number  of  predicted  pressures  within  i2  percent  of  their 
measured  values  increased  from  30  to  80  percent  compared  to  the  basic 
model . 

The  results  indicate  that  the  model  can  be  applied  to  inlet  pressure 
data  having  average  turbulence  levels  up  to  .04.  The  assumption  that  the 
pressures  have  a  normal  distribution  may  preclude  its  application  to 
those  situations  where  there  is  substantial  separation  and  planar  wave 
phenomena. 

The  relative  benefits  of  filtering  is  a  tunction  of  cut-off  filter 
frequency.  As  a  cut-off  filter  frequency  increases,  tfie  standard  devia¬ 
tion  of  the  difference  between  predicted  and  measured  pressures  increases 
and  the  number  of  pressures  agreeing  within  ±2  percent  of  their  measured 


values  decreases.  The  application  of  the  method,  including  map  averaging, 
may  be  limited  to  cut-off  filter  frequencies  of  less  than  1500  HZ. 

The  improved  method  has  been  shown  to  provide  predicted  pressure 
contour  maps  that  are  in  substantial  agreement  with  measured  contour  maps 
based  on  both  forty  and  eight  measured  turbulence  values. 
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